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We consider magnetotransport in high-mobility 2D electron gas, Oxx ^ 1, in a non-quantizing 
magnetic field. We employ a weakly chiral network model to test numerically the prediction of the 
scaling theory that the transition from an Anderson to a quantum Hall insulator takes place when the 
Drude value of the non-diagonal conductivity, CTxy, is equal to 1/2 (in the units of /h). The weaker 
is the magnetic field the harder it is to locate a delocalization transition using quantum simulations. 
The main idea of the present study is that the position of the transition does not change when a 
strong local inhomogeneity is introduced. Since the strong inhomogeneity suppresses interference, 
transport reduces to classical percolation. We show that the corresponding percolation problem is 
bond percolation over two sublattices coupled to each other by random bonds. Simulation of this 
percolation allows to access the domain of very weak magnetic fields. Simulation results confirm 
the criterion cTxy = 1/2 for values a^x ~ 10, where they agree with earlier quantum simulation 
results. However for larger ctxx we find that the transition boundary is described by Gxy ~ o'^^ with 
K ~ 0.5, i.e., the transition takes place at higher magnetic fields. The strong inhomogeneity limit 
of magnetotransport in the presence of a random magnetic field, pertinent to composite fermions, 
corresponds to a different percolation problem. In this limit we find for the delocalization transition 
boundary (Jxy ~ ctS-e'- 

PACS numbers: 72.15.Rn; 73.20.Fz; 73.43.-f 



I. INTRODUCTION 

Anderson localization is a single-particle phenomenon. 
Nevertheless, the scaling theory of localization^ which 
yields a profound prediction, full localization of all states 
in two dimensions, was formulated in terms of conduc- 
tivity of electron gas, a. Similarly, the extension^ of the 
2D scaling theory to a finite magnetic field is formulated 
in terms of components, Oxx and cr^y^ of the conductivity 
tensor of electron gas. Scaling equations describing the 
evolution of these components with the sample size, L, 
have the form 



1 



91ni 
da 



-aL2?e-2-^-cos(27ra,,), (1) 



aini 



= -ol^Ve--'-'^''- sin(2^a,,), 



(2) 



where 2? is a dimensionless constant. Drude values of a^x 
and a,j.y at size, L, of the order of mean free path, I, are 
given by 
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where ctq — kpl, kp is the Fermi wavevector, ujc is the 
cyclotron frequency, and r is the scattering time. These 
values serve as initial conditions to Eqs. ((!]) and ([2|). 
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termine the energies of delocalized states. 
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FIG. 1: (Color online) (a) Energy position of delocalized 
state, Eo, versus magnetic field, Wc, as predicted by Eq. Q. 
The curve, Eo{ujc), separates the phases with quantized Hall 
conductivities, axy ~ and a^y = 1. Cartoons illustrate elec- 
tron trajectories with restricted geometry in both phases; the 
edge state (red) is present in the upper cartoon and absent 
in the lower cartoon, (b) The predicted modification of the 
form of the flow diagram'' is illustrated schematically. 



The most nontrivial consequence of Eq. ^ is that it 
predicts levitation of delocalized states in weak magnetic 
fields ujcT <C 1, see Fig. [T^. In such fields it takes the 
form En — {n + ^)h/uJcT'^. In physical terms this means 
that a high-mobility electron gas with zero-field Drude 
conductivity axx = Epr/h 3> 1 exhibits a very strong 
sensitivity to a weak magnetic field. 
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as the temperature is decreased and quantum interfer- 
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ence effects become important. Phenomenon of levitation 
was predicted by Khmelnitskii^ even before Eqs. (H]) and 
^ were put forward (see also Ref. 0). Subsequently, it 
was also observed experimentally by several groups^^— . 

A. Physical interpretation of Eq. ([l]) 

The starting point in derivation of scaling equations 
^ and ([2]) was a cr-model with topological terni2.. It 
is desirable to understand physical processes underly- 
ing these equations. The first term Eq. ^ comes 
from Aharonov-Bohm phase action of the magnetic field. 
It describes that two paths corresponding to the same 
scatterers but different sequences of scattering events 
interfere even in the presence of the Aharonov-Bohm 
phases. The interpretation of the second term in Eq. 
([1]) is transparent in the limit of classically strong mag- 
netic field, > 1, where cos{2TTaxy) assumes the form 
cos{2TTEp/hujc), which is simply the field-induced modu- 
lation of the density of states. The origin of modulation 
is the emergence of Landau levels. On the other hand. 
Landau levels reflect the orbital action of the magnetic 
field, i.e., the fact that, with a certain probability, an 
electron can complete a Larmour circle with radius, i?^, 
without being scattered away by disorder. Thus the in- 
terpretation of the right-hand side of Eq. ([T]) is that the 
phase and orbital actions of the magnetic field compete 
with each other. 

Unlike strong fields, the interpretation of the cosine 
term in Eq. ([T} in weak fields, ^ 1, is much less 
transparent. In this limit we have (Jxy — ctqWcT in the ar- 
gument of cosine. The cosine term can be also rewritten 
as 



cos[2Tr{kpl){uJcT)] ~ 




where I b is the magnetic length, and <l>o is the flux quan- 
tum. For comparison, in the strong-field limit, the cosine 
term can be cast in the form 

cos(2^Bi?2/$o). (7) 

Comparing this expression to the last cosine in Eq. (jH]) 
suggests that in weak fields the role of the Larmour ra- 
dius is taken by the mean free path I. Note that I does 
not depend on magnetic field. Then the following ques- 
tion arises: what physics causes the orbital action of the 
magnetic field to manifest itself in the scaling equations 
in the weak-field limit? A possible way to unveil the 
orbital action is to adopt a cartoon picture where an 
electron moves not in a random potential but rather in 
a periodic background, say, on a quadratic lattice, as in 
seminal paper Ref. hj. Then we have to assume that 
the lattice constant, Z, is set by disorder. In this cartoon 
the orbital action will be encoded into the structure of 
the Bloch wavefunctions of electrons. It is the structure 
of the Bloch wavefunctions that leads to edge states in 



the presence of boundariesi^,. Note that the structure 
of the Bloch wavefunctions in a magnetic field depends 
crucially on the number of flux quanta through the unit 
cell, which, upon identifying I with the lattice constant, 
is the argument of the cosine in Eq. ([6]). Then the fac- 
tor, exp[— 27rA;jrZ], in front of the cosine in Eq. ([T]) has a 
meaning of degree to which a realistic random potential 
can be viewed as a periodic. Indeed, this factor can be 
interpreted as a probability for a realistic diffusive elec- 
tron to execute the same loop of length ^ I more than 
once. 

Obviously, realistic disordered system does not have 
any built-in spatial periodic structure. In view of the 
lack of a transparent interpretation of topological term 
in the weak-magnetic- field limit, it is important to check 
numerically whether or not some discrete value of the 
magnetic field of the order oi B ^ ^o/l^ causes delocal- 
ization transition and formation of edge state in a random 
potential. This was a subject of the papers Ref. [T3- In 
these papers, a network model describing a weakly chi- 
ral electron motion was introduced. The position of the 
quantum delocalization transition was established from 
the conventional transfer-matrix simulation of the trans- 
mission of the network. It was demonstrated that up 
to kpl ^ 10 the above estimate for the transition field 
applies. However quantum simulations become progres- 
sively complex in the limit of vanishing field. 



B. Delocalization transition with strong spatial 
inhomogeneity 

In the present paper we establish the position of de- 
localization transition indirectly. The underlying idea of 
our approach is that when a strong spatial inhomogene- 
ity is introduced into the quantum network, interference 
effects become progressively irrelevant, in the sense, that 
amplitudes of each two interfering paths typically differ 
strongly. Then the problem of the transport through a 
network reduces to the classical percolation. Most im- 
portantly, while the inhomogeneity-induced suppression 
of quantum interference leads to a strong reduction of the 
localization radius, the position of the transition remains 
unchanged. At the same time, classical simulations can 
be extended to much weaker magnetic field. The main 
outcome of our simulations is that for very weak fields 
(very small lOct) or high electron energies (very large kpl) 
the transition field is higher than ^>o/^^j namely 

B^^ikJT, (8) 

where k is close to 1/2. In terms of the flow diagram 
of the quantum Hall effect'^ the result Eq. ([8]) translates 
into the prediction that for Uxx > 10 the vertical flow 
line in Fig. [T)d deviates from axy = 1/2 to the right. 

In the present paper we also study levitation of delo- 
calized states in a vanishing average magnetic fleld, but 
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in the presence of a strongly fluctuating random mag- 
netic field. There is a notion that electron density varia- 
tions near the half-filling, v = 1/2, of the lowest Lan- 
dau level reduces to random magnetic field acting on 
composite fermions-l^''^- . It is also possible to realize an 
inhomogeneous magnetic field, acting on 2D electrons, 
artificiallyii"— . Different aspects of electron motion in 
random magnetic fields have been studied theoretically 
in Refs. lUll 

Fractional quantum Hall transitions can be associ- 
ated with quantization of cyclotron orbits of composite 
fermions. In this sense, fractional quantum Hall tran- 
sitions are the counterparts of delocalization transitions 
of electrons. Then the question arises: whether a de- 
localization transition for electrons at vanishing mag- 
netic fields has its counterpart for composite fermions 
at vanishing |;/ — 1/2|. At such filling factors, compos- 
ite fermion " feels" very weak average magnetic field. On 
the other hand, local fiuctuations of electron density give 
rise to a very strong random magnetic field, acting on 
composite fermion. For this situation we reduce the de- 
scription of magnetotransport to a different percolation 
problem. For critical values of filling factors we obtain 

\v-i/2\^[kAy 
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II. WEAKLY-CHIRAL NETWORK MODEL 
A. Description 

For completeness we remind the construction of the 
network model introduced in Refs. [l3 to describe quan- 
tum electron motion in a weak magnetic field. This con- 
struction is illustrated in Fig. [5] and consists of three 
steps: 

{i) We restrict electron motion by introducing forbid- 
den regions, Anm, (grey areas in Fig. [2]), which are not 
accessible for electron. Then electron moves in both di- 
rections along the links, which are the white regions, sep- 
arating Anm- The links join each other at the nodes, 
shown in Fig. [2] with brown full circles. 

(ii) We forbid forward and backward scattering at the 
nodes. This allows to parameterize the node scattering 
matrix, S'g, by a single parameter, q, as follows: 



(9) 



where Zi are the amplitudes of incoming and outgoing 
waves, see Fig. [21 

(m) We incorporate backscattering of electron moving 
along the link. The probability of backscattering is p, 
so that the corresponding scattering matrix, Sp, has the 
form 
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FIG. 2: (Color online) Left: Restricted electron motion over 
point contacts and bend-junctions is illustrated; An,m. are the 
centers of forbidden regions. Green line shows a minimal loop 
that can be traversed in clockwise and anticlockwise direc- 
tions. Right: Scattering matrices at the node and at the link. 



where Z\ and Z^ are amplitudes of incident waves, 
whereas Z\ and Z^ are amplitudes of refiected waves, 
see Fig. H 

To establish a correspondence with physical parame- 
ters, we identify the lattice constant with the mean free 
path, I. To relate the parameter q with the magnetic 
field, one can express, following Refs. |40| - |42| . the Hall 
resistance of the node, Rn, via the elements of matrix 
Sg, Eq. (HI): 

r + (1 - (lY 

In the absence of magnetic field Rh vanishes, indicating 
that (1/2 — q) is a measure of magnetic field, which is 
also the degree of preferential scattering to the left over 
scattering to the right. For a realistic electron moving 
a distance I in a magnetic field, this degree is WcT, thus 
allowing the following identification: 

i - g = WcT. (12) 

Note that electron motion over links and scattering at 
nodes in Fig. [2] models adequately the transport in the 
Boltzmann limit even at p = 0, i.e., without backscat- 
tering. However, a specifics of the diffusive motion with 
p = is that it does not allow quantum weak localiza- 
tion corrections. Indeed, weak localization corrections 
originate from the trajectories on the network for which 
an electron, starting from a certain link, returns to the 
same link with opposite direction of velocity (coherent 
backscattering). At p = the electron still can return 
to the same link, e.g., by encircling one forbidden region, 
but its velocity will be the same as the initial velocity. 
Finite p gives rise to weak localization. An example of 
an elementary loop providing coherent backscattering is 
shown in Fig. [21 The probability of this loop is 

r^p{l-pr[q{l-qf. (13) 

For realistic electron the return probability is {kpl)~^. 
This allows us to identify the parameter p as 
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FIG. 3: (Color online) Boltzmann transport on the p-q net- 
work. Rate equations Eq. (|15p relate the probabilities, 
Pi(m, n; t), to find electron on the corresponding half-link ad- 
jacent to the node with coordinates (m, n) and at time in- 
stances t and t + T. 



To conclude the construction of the quantum network, we 
assume as usually that random phases are accumulated 
in course of propagation along the links. This convention 
is non-trivial in the weak field limit. Indeed, as we dis- 
cussed in the Introduction, the delocalization transition 
is expected when the magnetic flux through a plaquette 
is of the order of flux quantum, $o- We will return to 
this point in Section IV D. 

The presence of two types of scattering processes, on 
the links and at the nodes, makes the Boltzmann de- 
scription of transport more complex. To develop this 
description, we turn to Fig. |31 It illustrates that the 
adequate variables to describe the Boltzmann transport 
are the probabilities, pi, i = I, .., 8, to find an electron on 
corresponding "half-link". In these variables, the closed 
set of rate equations reads 



plpeim - + pp2{Tn,n;t), 

q]P3 {m, n; t) + qprim, n; t), 
p]psim,n- l;t) + ppi{m,n;t), 
- q]p5{m,n;t) + qpi (m, n;t), 
p]p2im + l,n;t) + ppe{m,n;t), 
q]p7{m, n; t) -f qpz{m, n; t), 
p]p4{m,n + 1-t) +pps{m,n;t), 
q]pi (m, n; t) + qp^[m, n; t). 

(15) 



Performing Fourier transform in time and coordinate do- 
mains and taking the limit of small momenta, fc, and 
frequencies, a;^, we find a diffusive mode —iuJk — Dk^, 
where D is given by 
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FIG. 4: (Color online) Definition of the q-bonds. Scattering 
scenarios (a), (b), (c), and (d) correspond to the absence of 
both q-bonds, presence of both q-bonds, presence of one right- 
diagonal q-bond, and presence of one left-diagonal q-bond, 
respectively. 



q = 0, where the electron circulates around forbidden re- 
gions, clockwise and anti-clockwise, respectively. In these 
limits, for small p the diffusion coefficient is proportional 
to p. From Eq. (fT7|) we also see that D — >■ in the 
strong-scattering limit, p — > 1, as could be expected. In 
the limit of weak magnetic field, (1/2 — q) ^ 1, and high 
mobility, p <C 1 , we have 
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The fact that the magnetic-field correction to D is ^ 
{lOctY is generic for classical magnetotransport. On the 
other hand, the negative classical correction to D due 
to finite p is model-specific, since p was incorporated to 
capture interference effects. The meaning of the prefac- 
tor, Z^/r, which emerges from the system Eq. ((T5|) in the 
course of Fourier transform, is that the electron travels 
the distance of the mean free path, I, during scattering 
time, T. Correct, within a number, prefactor and mag- 
netic field dependencies of D indicates that the network 
model captures properly the magnetotransport in high- 
mobility electron gas in the Boltzmann limit. 

According to the scaling theory of localization, the 
knowledge of the Boltzmann transport coefficient should 
be sufficient to predict the position, Eq. (j4]) , of the quan- 
tum delocalization transitions, which in the limit of weak 
fields takes the form £^0 ^ filujcj^ . In the language of the 
network model this translates into the linear dependence. 
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p^-^q. 



(18) 



As discussed above, the diffusion coefficient is finite even 
at p = 0, except in the "strong-field" limits, g = 1 and 



Whether or not this prediction is valid can be estab- 
lished only by quantum numerical simulations. Espe- 
cially important is the limit, q — > 1/2, which corresponds 
to vanishing magnetic fields where a strong levitation is 
expected. Unfortunately this limit is the hardest to sim- 
ulate. This is because the localization radius to the left 
and to the right of the delocalization transition is huge. 
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ln(^//) = TT tJg ~ -K /p . This was a limitation of the 



quantum simulations reported in Refs. 
smallest value of p was p = 0.1. 



Il4, where the 



B. From quantum delocalization to classical 
percolation 

There is another, indirect, way to find the critical p — 
q boundary, bypassing quantum simulations, namely, to 
take the limit of strong disorder. By a limit of strong 
disorder we mean that local values pi and qi are strongly 
spread around averages p and q with distributions 
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Unlike the quantum case, where pi and qi were the same 
for all links and nodes, with distribution Eq. (fTO]) scat- 
terers on the links reflect fully in p percent of the cases, 
and transmits fully in the rest (1 — p) percent of cases. 
Similarly, according to Eq. (120]) . the nodes deflect only 
to the right in q"^ percent of the cases, deflect only to 
the left in (1 — q)^ percent of the cases; in the remaining 
2q{l — q) percent of the cases the deflection takes place 
both to the left and to the right depending on the in- 
coming channel, see Fig. |H The advantage of the strong 
disorder limit is that the quantum interference effects are 
irrelevant. The simplest way to see this is to turn to the 
elementary interference process illustrated in Fig. [2] If 
return to the origin is allowed for the clockwise direction, 
then it is forbidden for the ant i- clockwise direction since 
Qi(l — Qi) is zero in the strong-disorder limit. 

In the absence of interference the transport reduces to 
the classical bond percolation problem. The reduction is 
achieved by replacing scattering matrices Eqs. Q and 
(ITOl) by bonds according to the following rules: 

(i) The realization in which pi = 1 corresponds to 
quantum-mechanical reflection of incoming waves from 
all directions. In the language of percolation this con- 
flguration corresponds to a bond installed between the 
neighboring forbidden regions, and An+i^m, i-e., 

horizontal bond in Fig. [5] Below we will refer to this 
bond as a p-bond. For conflgurations with pi = the 
p-bond between the neighboring forbidden regions An^m 
and A„+i_m is absent. 

(m) The scattering matrix Sq is replaced by a pair of 
bonds (we refer to them as q-bonds), installed between 
the forbidden regions An,m and An±i,m±i, i.e., diagonal 
bonds in Fig. [21 Both q-bonds are absent. Fig. 0^, if 
the node deflects only to the left. Probability of this 
realization is Pa = {1 — q)^, as follows from Eq. ([20l) . 
Deflection only to the right corresponds to two crossed 
q-bonds present. Fig. This happens with probabil- 
ity Pb = q^ ■ The situation when right-diagonal q-bond 
is present while the left-diagonal q-bond is absent corre- 
sponds to the scattering scenario in Fig. 2};. The opposite 
scattering scenario. Fig. HJl, translates into left-diagonal 




FIG. 5: (Color online) Limit of strong disorder. The centers 
of forbidden regions, An^m and A„,m-i, are connected by the 
p- bond, while the centers of forbidden regions, An-\^m-i and 
An^rn, are connected by a q- bond. The delocalization tran- 
sition corresponds to the percolation threshold on the lattice 
consisting of p- and q- bonds. 



q-bond present and right-diagonal q-bond absent. The 
two latter bond conflgurations have equal probabilities, 
P, = Pd = q{l-q). 

Quantum-mechanical delocalization transition in the 
limit of strong disorder corresponds to percolation over 
p- and q-bonds, see Fig. [5] At the threshold of percola- 
tion p- and q-bonds form an inflnite cluster. At the same 
point, the waves propagating along the links in both di- 
rections and scattered at the links and at the nodes form 
an edge state. Threshold (9,^) values lie on a critical 
line of transitions on & q — p plain. Crucial for us is the 
relation between the points of this line and the positions 
of quantum delocalization transitions. In this regard it 
is important to relate the quantum matrix Sq to the ma- 
trices describing the different classical scenarios shown in 
Fig. m Setting q = we get 



Sa - 



-1 ON^ 
10 
1 

1000/ 



(21) 



which describes the scattering in Fig. |4^. Scattering 
scenario in Fig. |3]d is described by the matrix 




10 
0-10 
0-10 



1 



(22) 



which emerges upon setting q 
matrix 



1 in Eq. To get the 



Sc 



-1 
10 
0-100 
10 



(23) 



one has to set g = in the first and third columns, and 
q = 1 in the second and fourth columns. Similarly, the 
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matrix 

/ -1 \ 

Sd^ 1' (2^) 

V -10/ 

corresponding to Fig. |4ji emerges upon setting q = 1 in 
the first and third columns, and g = in the second and 
fourth columns. 

Matrices Sa and Sb provide non-zero Hall resistances, 
i?^ = — 1, i?^ — 1, while for Sc and Sd we have Rj^ — 
R% = 0, i.e., the Hall resistances are zero. The net Hall 
resistance is thus determined by 

PaR'h + PbRn ^ Pb - Pa, (25) 

which should be proportional to the magnetic field, 
(1/2 — q). The other relations between the probabili- 
ties of different scattering scenarios are normalization. 
Pa + Pb + Pc + Pd = 1, and obvious symmetry, Pc — 
Pd- These relations do not fix all probabilities uniquely. 
There is a profound physical reason for this ambiguity. 
Indeed, the net Hall resistance can be zero even if nodes 
locally deflect either to the left or to the right provided 
that Pa — Pb- This corresponds to the situation when a 
random magnetic field with zero average acts on electron, 
so that the time reversal symmetry is broken even in the 
absence of an external field. Such situation is generic for 
composite fermions, as was discussed in the Introduction. 
In addition to the probability assignment 

Pa^il-qf, Pb = q\ Pc = Pd = qil-q).. (26) 

dictated by Eq. ((20)) and described above, one can 
choose, e.g., 

Pa = l-q, Pb = q, Pc = Pd = 0, (27) 

when the electron scatters only to the left or only to the 
right from all incident channels. Obviously, for the latter 
assignment the magnitude of the random magnetic field 
is stronger than for assignment Eq. (|26|) . Finally, the 
physical situation when the time reversal symmetry is 
preserved in zero external magnetic field corresponds to 

Pa = 1-2(7, Pb = 0, Pc = Pd = q. (28) 

Three variants, Eqs. ([211), (113, and (US]), define three 
different percolation models, which we denote as A, B, 
and C, respectively. Results of numerical simulations of 
these models are reported in the next section. 

In conclusion of the present section we would like to 
draw a contrast between the classical limits of 4 x 4 
scattering matrix Sq and of 2 x 2 scattering matrix Sp, 
Eq. pn]) . Unlike the scattering matrix Sq, there is no am- 
biguity in taking the strong-disorder limit of the 2x2 link 
matrix Sp because this limit corresponds to the presence 
or absence of a single bond. In this regard, note that, in 
fully chiral network model by Chalker and Coddington'''^, 



scattering at the nodes is also described by a 2 x 2 scatter- 
ing matrix. The limit of strong disorder corresponds to 
the presence or absence of a single bond between the cen- 
ters of the squares A„_,„. Taking a strong-disorder limit 
in the Chalker-Coddington model reduces the quantum 
problem to conventional bond percolation on a square 
lattice. The position of the percolation threshold and 
the quantum delocalization transition certainly coincide, 
while the localization length in the strong-disorder limit 
is smaller—. 



III. SIMULATION PROCEDURE AND 
RESULTS 

In simulations performed, disorder realizations corre- 
spond to presence or absence of p- and q-bonds. In each 
realization, probabilities of p- and q-bonds are specified 
by the rules formulated above. Convention for the p- 
bonds, connecting counterpropagating links of 7i4-m-odd 
and n + m-even sublattices is the same for all three mod- 
els. Conventions for q-bonds are different for the models 
A, B, and C. These conventions are specified by Eqs. 
(PS)) . (1271), and (^51) . respectively. The main peculiar- 
ity of the simulations that complicates the trajectories 
stems from arrangement of pairs of q- bonds at the nodes. 
Namely, for different directions of approach to the given 
node the outcomes of passage are correlated. These cor- 
relations are illustrated in Fig. l5] The models A, B, and 
C differ by the weights with which different outcomes, a, 
b, c, or d, Fig. l5l are allowed. 

The size, L, of the samples used ranged between 500 
and 10000, where our unit of distance is half a link. For 
the largest system, we average over 10^ disorder real- 
izations. This number increases with decreasing size so 
that we keep a roughly constant CPU effort per size. 
To locate the position of the percolation threshold we 
searched for trajectories connecting two opposite faces 
of a square sample (periodic boundary conditions were 
imposed in the perpendicular direction). As in Ref. |45| . 
for a given realization, the two-terminal conductance be- 
tween the opposite open faces was identified with the 
number of such spanning trajectories. Different disorder 
realizations generate the conductivity distribution with 
average, <j{q,p, L). 



A. Phase diagrams 

To determine the critical boundary for each of the three 
models considered, PA{q), Peiq), and pc{q), we select a 
set of values of the turning probability q and then scan 
for many values of the probability p. For each size L, 
we represent the conductance as a function of p on a 
logarithmic scale and fit the points near the maximum 
of the conductance by a Gaussian. We then plot the 
position of the peak as a function of and extrapolate 
to infinite size. We found empirically that this fitting 
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procedure is of high quality for the three models. 

In Fig. |5]we represent the critical lines obtained for the 
three models: the upper curve corresponds to model B, 
the curve in the middle to model A, and the lower curve 
to model C. The straight line corresponds to p = 1/2 — q. 
The solid dots are the results of the quantum simulations 
in Ref. The inset shows the same critical lines as in 
the main panel at larger scale near the point (1/2, 0). 




0.0 0.1 0.2 0.3 



0.4 0.5 



FIG. 6: (Color online) Critical lines for the three model con- 
sidered: A (red, middle curve), B (blue, upper curve) and 
C (green, lower curve). The dots are the results of quantum 
simulations^''. Inset shows the same three critical curves at 
at higher scale near the point (1/2,0). 



We note that the curves tend to the point (1/2,0) as 
a power law with power higher than linear. To gain fur- 
ther insight, we analyze in detail the shape of the phase 
boundary at small probabilities p. In this regime qc is 
close to 1/2 and we expect a relation of the form 



p oc 



""2 



(29) 



In Fig. [7] we show p versus q — 1/2, on a double log- 
arithmic scale, for the three models A (middle set of 
points), B (upper set) and C (lower set). The straight 
lines are linear fits to the corresponding points. Their 
slopes are 7^1 = 1.994 ± 0.001, 7s = 2.464 ± 0.001 and 
7c = 2.286 ± 0.001. The errors quoted are the statisti- 
cal errors; systematic errors are also present, since it is 
impossible to include all finite-size effects. 




FIG. 7: (Color online) Critical lines for the three model con- 
sidered A (red, middle curve), B (blue, upper curve) and C 
(green, lower curve) on a double logarithmic scale near the 
point (1/2,0). 



B. Conventional percolation behavior away from 
q = 1/2. 

Quantum simulations in Ref. 14 demonstrated that 
the delocalization transition along the boundary, p{q), 
belongs to the quantum Hall universality class. In par- 
ticular, it was demonstrated that for the first three black 
dots in Fig. [6] which correspond to q < 0.3, the critical 
exponent is close to 7/3. Introducing strong disorder sup- 
presses the quantum interference. We expect that quan- 
tum percolation at a given {q,p) reduces to classical per- 
colation and the critical exponent = 7/3 is replaced by 
its classical value = 4/3. For the Chalker-Coddington 
model, in which the position of delocalization is fixed 
at average value of the disorder potential, the crossover 
from 7/3 to 4/3 was tested in Ref. ^45. In this section we 
demonstrate that, away from the point p = 0, g = 1/2, 
the p{q) boundary established above indeed corresponds 
to the divergence of the localization length with exponent 
1/ = 4/3. 

The determination of the critical exponent is based 
on the fact that near {q,p) — {qc,Pc) the conductance 
is a function of a single argument, {p — pc)L^^^ (ver- 
tical scan) or {q — qc)L^^^ (horizontal scan). Exactly 
at (qcPc) the conductivity assumes the universal value, 
(70 = 0.361404 . . . found by Cardy^^,, which should be the 
same for the entire boundary. 

The scaling analysis was performed for all three mod- 
els. In Fig. [5] we present results for the model A a.t a 
particular critical point (0.3,0.23392). Overall, the scal- 
ing confirms that v = 4/3 both for vertical and horizontal 
scans. Particular feature about the scaling data is that 
the widths of scaling functions are slightly different for 
the vertical and horizontal scans. We have also found 
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that there is a small size effect precisely at the boundary 
which is well described by the expression 



where a is a constant. 
0.4 



L3/4' 



0.3 

■> 

t3 0.2 

c 
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" 0.1 



0.0 



(30) 
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FIG. 9: (Color online) Conductivity as a function of p for 
model A along the line q = 1/2. The lateral sample sizes are: 
500 (cyan diamonds), 1000 (blue down triangles), 2000 (green 
up triangles) , 4000 (red circles) and 8000 (black squares) . In- 
set shows the position of the peaks as a function of i on a 
double logarithmic scale. 



FIG. 8: (Color online) Scaled conductivity as a function of 
the probability difference with the critical point multiplied 
by L^'^*. The solid symbols correspond to vertical scans and 
the empty symbols to horizontal scans crossing the critical 
point (0.3, 0.23392) of model A. The lateral sample sizes are: 
500 (cyan diamonds), 1000 (blue down triangles), 2000 (green 
up triangles) and 4000 (red circles). 



C. Behavior of the localization length at zero field: 
models A and C. 

We now turn to the behavior of the localization length, 

at zero magnetic field, q = 1/2. For each model, the 
dependence ^{p) at small p is determined by a peculiar 
behavior of the corresponding delocalization boundary 
established in subsection A. It is also very important that 
^(p) is equally affected by the second, complementary, 
boundary in the domain q > 1/2, which is the mirror 
image of the boundary in Fig. [51 

For moderate p, when the boundaries were straight 
lineai^, the presence of the second boundary leads to the 
enhancement of ^. On the contrary, we will see that in 
the limit p — 0, the fact that both boundaries for a given 
model are almost horizontal leads to a shortening of ^. 

In Fig. ini we plot the conductivity for model A as 
a function of p for q = 1/2 and for several values of 
the system size: 500 (diamonds), 1000 (down triangles), 
2000 (up triangles), 4000 (circles) and 8000 (squares). It 
is seen that the curves a{p) for different system sizes have 
similar shapes and are even-spaced along the logarithmic 
horizontal axis. Thus we expect scaling and behavior ^ ^ 
p^'^^ as a consequence. A practical procedure to infer 
from the data in Fig. [5] is based on the dependence, pl 



versus L, where is the position of the maximum of 
the conductivity for a given L. In the inset of Fig. [HI we 
plot pl{L) in a double logarithmic scale. We see that pL 
follows the dependence 

PL=bL-P, (31) 

where h and /3 are model-dependent constants. For model 
A, we found = 1.51±0.02 and 6 = 4.1±0.8. Eq. jSl]) 
and the fact that is very close to 3/2 suggest that, to 
achieve scaling, the data in Fig. [HI should be replotted 
versus pl?!"^ . The result of this replotting is shown in Fig. 
[TOl It is seen that the overlap is excellent, yielding the 
critical exponent — 2/3. This should be contrasted to 
the behavior ^(p) cx (p — Pc)^^^"^ at an?/ non-zero pc- We 
conclude that at = the divergence of the localization 
length with p is much slower. 

Finally, for the model C the plots ct(p) as a function of 
p do not exhibit maxima. As shown in Fig. Ill[ where we 
plot ct(p) versus pl? ^ a very good overlap is achieved 
for vc = 4/7. 

In conclusion of this subsection we note that the con- 
ductivity for all three models tends to 2cto — 0.722808... 
as p — > 0. The reason is that the value cr = ctq at the 
threshold is the property of a single critical point^^ . By 
contrast, in our case two critical lines merge at the point, 
q = l/2,p = 0. 



D. Behavior of the localization length at zero field: 
model B. 

Scaling analysis of the data for model B reveals slightly 
different behaviors for the domains of "moderate" p > 
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FIG. 10: (Color online) Conductivity as a function of pL^^^ 
on a logarithmic scale for model A along the line q = 1/2. 
The sample sizes are: 500 (cyan), 1000 (blue), 2000 (green), 
4000 (red) and 8000 (black). 



FIG. 12: (Color online) Conductivity as a function of pL7^'^ 
along the line q = 1/2 on a logarithmic scale for model B. 
The sample sizes are: 1000 (blue), 2000 (green), 4000 (red), 
6000 (cyan), 8000 (black) and 10000 (magenta). Inset: the 
same as main plot with low-p data included; conductivity is 
plotted versus pL^^^. 




100 



FIG. 11: (Color onhne) Conductivity as a function of pL'^'^ 
on a logarithmic scale for model C along the line q = 1/2. 
The sample sizes are: 500 (cyan), 1000 (blue), 2000 (green), 
4000 (red) and 8000 (black). 



10""* and "truly critical" p < 10~^. For the first domain, 
from the position of peaks we find /Sg = 1.77 ± 0.04. 
This suggests that t'e = l//?e~4/7. Note however that 
replotting the conductivity versus pL^^^, see Fig. [T^l 
does not lead to overlap as good as for the model A. 
Moreover, in the second domain p < 10~* a good overlap 
is achieved for the exponent 2/3, i.e., the same as in the 
model A. This is illustrated in the inset of Fig. [T^l This 
indicates that for the model B the true critical region 
is quite narrow. Such a delicate behavior of ^(p) for the 
model B might indicate that the critical boundary p{q) in 
this model also changes the behavior in the truly critical 



IV. DISCUSSION 



Position of boundaries 



region p 



< 10- 



It is seen from Fig. [6] that the boundaries, pa{q) and 
Pc{q), almost coincide in the entire domain, < q < 1/2. 
Overall, these boundaries are in a gree ment with the re- 
sults of quantum simulation Ref. [ij shown with black 
dots. It is also seen that the boundary, psil), goes sig- 
nificantly higher. In particular, at q = 0.25, pe exceeds 
PC almost twice. On the physical level, this means that, 
for a given average magnetic field, the formation of edge 
states requires a longer zero-field mean free path, I, for 
model C than for model B. In other words, formation 
of edge states happens easier when a random magnetic 
field is present. To gain a physical insight why this is so, 
consider electron motion in random magnetic field. Lo- 
cal value of the field changes its sign in space, while the 
average field, (1/2 — g), is much smaller than the absolute 
value of the local field. Then electron trajectories are ei- 
ther circles inside the regions where the field maintains 
its sign, or snake states, propagating along the bound- 
aries of these regions, i.e., along the contours with zero 
local field. Then it is apparent that a weak disorder does 
not affect this picture. If, on the other hand, the mag- 
netic field, (1/2 — q), is uniform, electron trajectories are 
big circles. Then a weak disorder will have a strong effect 
by deflecting electron before it completes a circle. The 
above two situations correspond to the models B and C, 
respectively, and explain why psiq) > Pcil)- In model 
A, random component of magnetic field is present, but is 
weaker than in model B. In this regard, the fact that the 
boundary PA{q) lies between psiq) and pc{q), also finds 
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its explanation. 

Although quantum simulations are impossible in the 
domain p ~ 0.01 five "quantum" data points in Fig. [6] 
follow pc{q) within the accuracy of quantum simulations. 
Then the confirmation of the scaling theory Eqs. (H]), ([2]) 
would be the linearity of the percolation boundary at 
q 1/2, see Eq. 

The most important outcome of the present simula- 
tion is the inset in Fig. |6l It is seen that at really 
small p ^ 0.01 and q close to 1/2 the behaviors of all 
three boundaries changes dramatically compared to their 
"bodies", namely, they become almost horizontal. All 
three boundaries have the form p ^ (1/2 — g)'' with 
7 > 2. This is in stark contrast to the prediction of 
scaling theory Eq. (|18p. which corresponds to 7 = 1. In 
other words, percolation results suggest that instead of 
the condition a^y — 1/2, the delocalization boundary is 
described by 

i-i 

CTxy ^ CTxx '' ■ (32) 

The latter condition can be also cast in the form Eq. 
([8]) with K = 1 — 1/7. We note that the crossover from 
'^xy — 1/2 to Eq. takes place at large Or^,^ ~ 10. 

In terms of the flow diagram of the quantum Hall effect^ 
this means that the upper part of the vertical flow line is 
bent to the right, as it is illustrated in Fig. ^p. 



B. Semi-analytical consideration 

To specify the distinct behavior of percolation bound- 
aries in vanishing average magnetic field they are plotted 
in Fig. [7] in the log-log scale. From the slopes we deduce 
the values, 7^ = 1.994, 70 = 2.464, and 7c = 2.286. To 
get a feeling why all 7- values are close to 2, below we 
present some semi-analytical arguments. We first turn 
to Fig. [5] and set p = 0. Then the lattice breaks into two 
quadratic sublattices with n ~\- m even and n + m odd, 
which are completely disconnected. None of them perco- 
lates if q, the percentage of bonds present in each sublat- 
tice, is less than 1/2. Finite p = Pc{q) allows percolation 
for g < 1/2 since p-bonds couple clusters from different 
sublattices. It is apparent that coupling of clusters by p- 
bonds is relevant if the typical distance, 1 / between 
two p-bonds become smaller than the localization length, 
i{q) = (l/2-g)-*/3. This yield a constrain that 7 < 8/3. 
This constrain is insensitive to the mutual correlations of 
q-bonds on the two sublattices. In fact, this correlation is 
absent in model A. Indeed, as follows from Eq. at 
q = 1/2 for model A we have Pa = Pb = Pc = Pd = 1/4. 
By contrast, for model B the probabilities at q — 1/2 
are P^ = Pb = 1/2, = = 0. This suggests that 
q-bonds on two sublattices are strongly (and positively) 
correlated. Namely, if there is a q-bond connecting two 
n -\- m even plaquettes at a given node, then there must 
be a q-bond connecting n + m odd plaquettes at the same 
node. On the other hand, the correlation of q-bonds at 




FIG. 13: (Color online) Vicinity of the point p = Q, q = 
1/2, of the phase diagram Fig. [T] A q-cluster of n -I- m odd 
sublattice (upper) and a q-cluster of n -I- m even sublattice 
(lower) overlap. A joint trajectory (thin blue line) is formed 
upon installing of a single p-bond. The blowup illustrates 
hybridization of trajectories on a microscopic level. 

a node in model C is negative: presence of one q-bond 
excludes the presence of the other. This is apparent from 
Fig. [2 

As p-bonds are switched on, clusters include sites from 
both sublattices, see Fig. [T31 In Ref. |I3 where the 
model A was considered, it was argued that 7=1. The 
argument was based on the following picture of the clus- 
ter growth upon increasing q: it was assumed that criti- 
cal q-clusters on a given sublattice grow by getting con- 
nected via additional q-bonds. Since the growth of clus- 
ters due to p-bonds takes place by connecting critical 
clusters from different sublattices, it was concluded that 
p- and q-bonds play equal roles in the growth of clusters, 
which immediately leads to 7 = 1. Present simulations 
suggest that in the close proximity of g = 1/2 this pic- 
ture fails, and the role of p- and q-bonds in approaching 
the percolation threshold is completely different, namely, 
the growth due to p-bonds is more efficient. This growth 
proceeds by p-bonds connecting the hulls of critical clus- 
ters from two sublattices, as illustrated in Fig. [T31 For a 
given g, the length of the hull is 

C[q)=e'\q)=[^-;^-qj ■ (33) 

For the model A, to achieve a percolation by adding p- 
bonds one should take into account that the hulls on two 
sublattices are uncorrelated. Then, a p-bond with one 
end on a hull from even sublattice will have the other end 
on the hull from odd sublattice with probability Cj^ . As 
a result, the percolation condition reads 

{Cie){vC) = 1, (34) 

where the second factor is the probability that there is at 
least one p-bond with one end on the critical hull from, 
say, odd sublattice. Eq. yields 7,4 = 2, which coin- 
cides with the simulation result. 
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Correlation of q-bonds at the nodes for models B and 
C leads to conclusion that corresponding critical clusters 
in two sublattices are also correlated. One consequence 
of this correlation is that it takes less p-bonds than in 
model A to connect critical hulls from two sublattices. 
As a result, 7e, 7c > 7^ = 2. On the other hand, the 
picture of percolation by connecting the critical clusters 
imposes the upper boundary 7 < 7/3 for both models B 
and C. This follows from the condition that there should 
be at least one p-bond per critical hull, i.e., pC> 1. Note 
that the constraint 7 < 7/3 is stricter than the constraint 
7 < 8/3, established above. 

Beyond the estimate, 2 < 7^, 7c < 7/3 we cannot 
come up with more accurate analytical values for these 
indices. We are not even able to establish which of them 
is bigger. This is because strong correlation between the 
hulls in both models B and C simplifies connectivity upon 
switching on p-bonds. On the other hand, this correla- 
tion prevents the expansion of the resulting cluster. 

As it was established in the previous section, in the 
domain of p < 10"'', behavior of ^{p) a,t q = 1/2 in the 
model B exhibits crossover from the critical exponent, 
4/7, to 2/3. To relate this peculiar behavior with the 
shape of the percolation boundary psiq) = (1/2 — q)^'^ , 
we invoke the argument of Ref. 47 which, in application 
to the p-q model, goes as follows. If the divergence of 
^ at the point, q — 1/2, p — 0, is characterized by ^ ^ 
(1/2 — q)^'^' along the q- direction and ^ ^ p^'^p along 
the p- direction, then the shape of the critical boundary 
is p ^ (1/2 — qf'^l"'^ , i.e., 7 = Vq/vp. Following this 
argument, crossover in the model B from Vp — vb — '^H 
to vti = 2/3 suggests that 70 and 771 merge in truly 
critical region. 

Overall, our numerical results suggest that in the truly 
critical domain, where 7.A ~ 78 ~ 2 and 7c ~ 7/3, the 
divergence of ^(p) for all three models is well described 
by the relation 

^^p-4/(37), (35) 

With regard to the argument of Ref. [43 this means that, 
for all three models, the exponent Vq is equal to 4/3, i.e., 
the same as for q away from 1/2. 

C. Relation to Ref. gl 

In Ref. IZsl . spin quantum Hall effect in bilayer and tri- 
layer systems was studied numerically, in order to trace 
the emergence of macroscopic metallic phase upon adding 
the third dimensior^^. The authors made use of the 
fact that in a strictly 2D system there is a mapping 
between the spin quantum Hall transition and classi- 
cal bond percolation^Si^. For bilayer systems the corre- 
sponding classical percolation is bond percolation on each 
layer (bonds connect the centers of plaquettes) , comple- 
mented with the possibility to switch layers with a proba- 



bility, pi , while passing each side of each plaquette. Phys- 
ically, in spin quantum Hall effect, an electron travels on 
each layer of the network in the same direction. In our 
consideration of the weak-field quantum Hall effect, an 
electron stays within a plane but each link of the square 
lattice represents two counterpropagating channels. For 
this reason there is a mapping between the simulation in 
Ref. and treatment of the model A in the present pa- 
per. Namely, pi in Ref. |^ should be identified with 
the backscattering probability p in the present paper, 
while the probability that the given bond is present in 
Ref. |45| . p, should be identified with our parameter q. 
Due to this mapping, critical behavior, pi{p) for small pi 
in Ref. H^, is the same as the behavior of critical line, 
p fx. (1/2 — q)^^, for small p in our model A. Also, the 
above semi- analytical calculation of 7^ is the same as 
proposed in Ref. H^. However, it should be noted that 
mapping between the network of Ref. |4^ and model A 
applies only for small »i . For larger pi the position of 
the boundary in Ref. |45| differs dramatically from that 
of the model A. 



D. Phases with higher a^y 

In fact, scaling theory predicts that all Landau levels, 
n, in Eq. (|4]), eventually levitate to £"„ -> 00, as the 
magnetic field is lowered. Our simulations do not cap- 
ture low-field transitions for n > 1. This is because we 
restricted our consideration to network with one channel 
per link. Within this description we were able to capture 
the Drude conductivity tensor of electron gas in a weak 
field, and weak localization effects. On the other hand, 
this description does not allow, in principle, to capture 
the phases with quantized Uxy higher than 1. 

As a final remark, we can underscore the difference 
of the scaling theory and our results as follows. The 
scaling theory predicts that electron gas experiences a 
delocalization transition in a magnetic field at which fiux 
into the area P is of the order of the flux quantum, $0- 
We find that, in the limit of kpl > 10, the interplay of 
orbital and phase actions of magnetic field causes the 
transition when this flux is much larger than $o- 
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